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The critical behavior of the random-field Ising model has been a puzzle for a long time. Different 
. theoretical methods predict that the critical exponents of the random-field ferromagnet in D dimen- 

sions are the same as in the pure (D ~ 2)-dimensional ferromagnet with the same number of the 
magnetization components. This result contradicts the experiments and simulations. We calculate 
£N| ■ the critical exponents of the random-field 0(A) model with the (4 + e)-expansion and obtain values 

' different from the critical exponents of the pure ferromagnet in 2 + e dimensions. In contrast to 

O ; the previous approaches we take into account an infinite set of relevant operators emerging in the 

problem. We demonstrate how these previously missed relevant operators lead to the breakdown of 
the (6 — e)-expansion for the random-field Ising model. 

The role of quenched disorder in condensed matter depends on its strength but even weak disorder can strongly 
^ • modify the nature of the phase transition. This effect is most prominent in the case of random- field (RF) disorder which 
breaks both the translational symmetry and the symmetry with respect to transformations of the order parameter. 
In particular, in some systems the arbitrarily weak disorder of this type destroys long-range order 0. The strong 
effect of weak disorder makes it difficult to apply the standard perturbative methods of phase transition theory to RF 
problems. Besides, it is much more difficult to solve exactly a random model than a pure one. As a result, a theory 
*h ■ of phase transitions in the presence of random fields is still absent. 
'■^J | The large number of systems with RF disorder provides a strong motivation to develop such theory. Some of these 
' systems are known for a long time. Examples are disordered antiferromagnets in the external mag netic field [§, 
binary liquids in random porous media || and vortices in disordered superconductors Q| . Recently a lot of attention 
i . was devoted to disordered liquid crystals ||] and liquid He-3 in aerogels [|| . In contrast to the pure systems even the 
question of ordering at low temperatures in the presence of random fields is nontrivial. The more difficult problem of 
t-H the critical behavior is still open in spite of two decades of investigations. 

After a lot of controversy the lower critical dimension of the RF Ising model was found |l| but the structure of 
the phase diagram and the nature of the transition to the ferromagnetic state are still unclear. Different theoretical 
approaches to the paramagnet-ferromagnet transition predict the dimensional reduction the critical exponents 
of the RF 0(A) model in D dimensions should be the same as in the pure 0(A) model in D — 2 dimensions. This 
prediction relates the three-dimensional RF Ising model to the one-dimensional pure system in contradiction to the 
presence of long-range order in the former and its absence in the later. Moreover, the high-temperature expansion 
shows that the dimensional reduction rule is invalid in any dimension less than the upper critical dimension 6. 
The most elegant derivation of the dimensional reduction || provides an (incorrect) exact solution of the RF 0(A) 
model at zero temperature. The failure of that approach was explained by Parisi on the basis of the complicated 
energy landscape of the model However, such explanation is insufficient for the renormalization group (RG) and 
1/A-expansion. It was conjectured |9| that RG fails because of non-perturbative corrections to the expansion in e. 
However, such corrections are not found and recent numerical results |Q suggest that the existing predictions for the 
q ' coefficients of the series in e are wrong jlO| . 

, A possible reason why the perturbation theory does not provide satisfactory results is the appearance of some 
additional relevant operators missed by the existing approaches. This phenomenon is responsible for the failure of 
the standard RG theory to predict the order of the phase transition in some systems . The possibility of a similar 
phenomenon in the RF 0(A) model was first suggested in Ref. [|l2| but no new values of the critical exponents were 
found. Recently a similar idea was used in Ref. pi but it also did not allow to calculate the exponents. Besides, the 
approach |r| involves delicate manipulations with the replica limit and its validity is not clear fl4|| . In the present 
paper we demonstrate that there are indeed some additional relevant operators in the problem although not those 
found in Ref. We calculate the critical exponents of the RF 0(A) model in 4 + e dimensions and show that 

they do not satisfy the dimensional reduction. We also demonstrate how the additional relevant operators lead to 
the breakdown of RG in 6 — e dimensions. Some phenomenological approaches |l5| ] allowed one to obtain the critical 
exponents different from the dimensional reduction prediction. However, their results based on different unproven 
assumptions contradict each other. An important breakthrough was made by Mezard and Young |l6| ] who considered 
the possibility of the replica symmetry breaking in the RF 0(A) model at large A. Unfortunately, the approach |l6} | 
did not allow the calculation of the exponents. The present paper contains the first systematic method that explains 
the failure of RG and allows us to find the critical exponents of the RF systems. 
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In the critical point the RF O(N) model can be described by the connected and disconnected correlation functions 
0: 

G con (q) = <[n(q)n(-q)] - [n(q)][n(-q)]> ~ <T 2+ "; G dis (q) - <[n(q)n(-q)]} ~ g"*^, (1) 

where n(q) is the Fourier component of the magnetization, the square brackets denote the thermal average and the 
angular brackets denote the disorder average. To calculate the critical exponents r\ and fj we develop the e-expansion 
near the lower critical dimension 4 following the line of Ref. |Tij] . Our starting point is the Hamiltonian of the RF 
O(iV) model 

H= I ^^a;[J^^n(x)^n(x)-^(h fc (x)n(x)) fe ], (2) 

J ft k 

where the unit vector of the magnetization n(x) has N components and are random fields with zero average, 
(hk 1 a{x)hq l /3{y)} = Hk5 a f38kqS(x — y). The Hamiltonian includes an infinite set of random anisotropies of different 
ranks. These contributions are allowed by symmetry and turn out to be relevant in the RG sense. The replica 
Hamiltonian has the form 

^=^€^^^^ + •■•1 (3) 

a fj, ab 

where a, b are replica indices, R(z) is some function, T the temperature and the dots denote the irrelevant terms. Near 
the zero-temperature fixed point the whole function R(z) is relevant. Indeed, to ensure the fixed length condition 
n 2 = 1 at each RG step, we ascribe the scaling dimension to the magnetization n. The scaling dimension of the 
temperature is —2 + O(e). The relevance of any operator is determined by the number of the derivatives in it and the 
power in which it contains the temperature. This shows that all operators Rk — J2ab( n a n b) k /T 2 are relevant in the 
same space dimensions. The functional RG analysis in 4 + e dimensions repeats the same steps as the (4 — e)-expansion 
developed in Ref. jl7) and results in the following RG equation for the disorder strength 

= = -eR{<t>) + {R"m 2 2R"(<f>)R"(0) - 

(N - 2)[4R(<j>)R"(0) + 2ctg<W)i?"(0) - (^] } + 0(R 3 ,T), (4) 

\ sm0 / 

where we define cos <f> = n a n;, to make the equation more compact. The RG equation is valid in the zero-temperature 
fixed point that describes the phase transition in the RF ferromagnet. The critical exponents ([!]) can be expressed 
via the RG charge R(<j>): 

r] = -2R"(0); fj = -2{N - 1)R"(0) - e. (5) 

Note that R"(0) < since G CO n(r) ~ r 2 ( A '~ 1 ) fl (°) must be limited. The solution of Eq. (^) can be found numerically 
with shooting. Equation (Q) has only a discrete set of the solutions which are everywhere finite. The region of possible 
initial conditions R"(0) is limited by the restrictions 

e/2(N - 3) > -R"(0) = r//2 > e/2(N - 2). (6) 

The first inequality follows from the Schwartz-Soffer inequality [jl8| and is the stability condition for the fixed point 
p7[. The second one can be derived from the third term of the expansion of the function R(<fi) at small <j>: 



{N-l)R"{0f R "(°),2 + R"(Q)e + 2(N-2)R"(0r ^^ 
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Since coefficients of Eq. (||) are singular at 4> = 0, tt, one has to use the expansions of R{4>) in powers of 
near those points at solving numerically the RG equation. Note that the function R(<j>) is non-analytical at <f> = 0. We 
have solved the RG equation at N < 5. It turns out that for 5 > N > 3 Eq. ([I|) has only one solution compatible with 
the condition (||). The critical exponents calculated for those N are given in Table I. One can see that the critical 
exponents are different from the dimensional reduction values 77(4 + e) = 77(4 + e) = rj purc (2 + e) = e/ (N — 2), where 
?7 P urc is taken from Ref. p9| . 
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The second of the inequalities (||) shows that for the XY model (N — 2) there is no appropriate solution of the RG 
equation. This can be explained by the fact Q that the RF XY model possesses quasi-long-range order in dimension 
3. Hence, it has a phase transition in 3 dimensions and thus its lower critical dimension near which the e-expansion 
should be developed is less than 4. 

It was claimed p0| that in the first order in e the critical exponents of the RF ferromagnet in 4 + e dimensions are 
the same as in the pure 0(A) model in 2 + e dimensions. That prediction could be obtained from Eq. (Q) with the 
wrong assumption that there is an analytical solution R{<j)) controlling the critical point. 

We have demonstrated the existence of an infinite set of relevant operators near 4 dimensions in the RF 0(A) 
model. A possible explanation of the failure of the (6 — e)-expansion is thus the appearance of additional relevant 
operators below some dimension D < 6. However, we show that another scenario takes place: an infinite set of relevant 
operators emerges in any dimension less than the upper critical dimension 6. Since the upper critical dimension is 
the same for the Ising and 0(A) models our approach allows us to consider both cases in the same way. 

Near 6 dimensions it is convenient to use the Ginzburg-Landau model with the random field h(x) 

H = J d D x[{Vm) 2 + 5 (m 2 ) 2 - h(x)m] (8) 
that can be described by the replica Hamiltonian 

H R = /AElVm^+^Kf-A^mA], (9) 

a a ab 

where m a are the replicas of the A-component magnetization. The standard power counting suggests that all operators 
which are relevant in 6 — e dimensions are included into Eq. (|9|). However, if it were so the dimensional reduction 
would be correct. The hint why the power counting fails is given by the theory of the metal-dielectric transition. It 
was argued that in the nonlinear sigma-model for the metal-dielectric transition one should include an infinite set of 
relevant operators J2l]] missed by the power counting. For any eigen-operator A k of the RG transformation the RG 
equation has the structure 

A A, 

' " v (C* + C*e s + o(e s ))A k + 0(A 2 ). (10) 
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The power counting is based on the sign of the constant Co and predicts that the operator A k is irrelevant as Cg < 0. 
However, if the ratio C\ /Cq grows up to infinity at large k then for any fixed e there is such k that the correction 
(jk e s j g g rea t er than Cq. Hence, the role of the operators A k depends on the signs of C\. In the 4 -theory without 
the random field this sign is negative for any operator which is expected to be irrelevant from the power counting 
p2| . This agrees with the success of the RG approach in that problem. On the other hand, we shall see that the signs 
of these constants are positive for some operators in the RF problem. This signals that additional relevant operators 
emerge. 

For the pure 4 -model one can define the rank of an operator with q derivatives and the pth power of the order 
parameter as r = p + q — 4. Since only the operators of the same rank can mix pl| , for the calculation of the anomalous 
dimension C^ e s of an operator with rank r k one can ignore the diagrams which produce the operators of the higher 
ranks. It is easy to check that the definition of the rank has to be modified in the RF problem to ensure that the 
operators of different ranks do not mix: r=p + q + 2t — 6, where t is the number of the different replica indices in 
the operator. 

To avoid the difficult problem of operator mixing we guess a set of relevant operators which do not mix with the 
other operators of the same rank up to the second order in e. It turns out that the following operators are relevant 
and do not mix with the other: A k = J2 a b[{m a — mj,) 2 ]'. At k > 1 their canonical dimensions 

d k c = (4- e )fc-6 + e (11) 

are positive and proportional to k. We shall see that the anomalous dimensions are negative and proportional to (ke) 2 . 
If one imposes the fixed length restriction m 2 = 1 so that the Ginzburg-Landau model reduces to the 0(A) model 
then the operators A k can be represented as the linear combinations of the random anisotropics B k = 6 (m a m;,) . 
This is natural since the operators B k are relevant in 4 + e dimensions as one could see above. 

To the first order in the quartic vertex g the anomalous dimensions of the operators A k are zero. To prove this 
we demonstrate that all diagrams with one vertex gm^ and one vertex Af? = [(m a — mj,) 2 ]* 1 either are equal to zero 
or produce operators of higher ranks. Indeed, any line of the Feynman diagram corresponds to the product of a 
momentum-dependent factor and two operators which differentiate the vertices at the ends of the line with respect 
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to replicas m a of the magnetization. For example, line 2 in Fig. 1 acts on the vertex as the differential operator 
82 = d/dm a . Line 1 in which the vertex A is inserted acts on the vertex A% b as the operator d\ = d/dm a + d/dmi,. 
Obviously, diA% b = 0. Hence, any diagram including line 1 is equal to zero. On the other hand, any diagram with 
one vertex g, one vertex A/., and without lines in which A is inserted produces an operator of a higher rank. Fig. 2 
shows the only nonzero diagram of the order g 2 that does not increase the rank of the operators Ak . Calculating this 
diagram one obtains the anomalous dimension 

= k(N + 2)e 2 /2(N + 8) 2 - e 2 [Nk(2k + 1) + 16fc 2 - 10fc]/2(7V + 8) 2 (12) 

The anomalous dimension (^) is negative and proportional to k 2 . Thus, at any fixed e one expects that the 
operators Ak with 1/e 2 < k are relevant but missed by the existing theoretical methods. Although the rigorous 
analysis requires consideration of all orders in e, the appearance of an infinite set of relevant operators in 6 — e 
dimensions is plausible since such set exists near 4 dimensions and since the alternative explanation of the failure 
of RG due to non-perturbative corrections is hardly compatible with the existing numerical results 

In conclusion, we have calculated critical exponents of the O(N) model in 4 + e dimensions and demonstrated that 
they do not obey the dimensional reduction. The failure of the dimensional reduction is related to the appearance of 
an infinite set of relevant operators. 
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the critical reading of the manuscript. This work was supported by the Koshland fellowship and RFBR grant No. 
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TABLE I. Critical exponents of the RF O(N) model. r\ and r\ are the exponents of the connected and disconnected correlation 
functions respectively. 



N 


3 


4 


5 


V 


5.5e 


0.79e 


0.42e 


n 


lO.le 


1.4e 


0.70e 



line 1 




line 2 



FIG. 1. A first order diagram contributing to the anomalous dimension. 
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FIG. 2. A second order diagram contributing to the anomalous dimension. 
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